Kinetic enhancement of Raman backscatter, and electron acoustic Thomson scatter 
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1-D Eulerian Vlasov-Maxwell simulations are presented which show kinetic enhancement of stim- 
ulated Raman backscatter (SRBS) due to electron trapping in regimes of heavy linear Landau 
damping. The conventional Raman Langmuir wave is transformed into a set of beam acoustic 
modes [L. Yin et al, Phys. Rev. E 73, 025401 (2006)]. For the first time, a low phase velocity elec- 
tron acoustic wave (EAW) is seen developing from the self-consistent Raman physics. Backscatter 
of the pump laser off the EAW fluctuations is reported and referred to as electron acoustic Thom- 
son scatter. This light is similar in wavelength to, although much lower in amplitude than, the 
reflected light between the pump and SRBS wavelengths observed in single hot spot experiments, 
and previously interpreted as stimulated electron acoustic scatter [D. S. Montgomery et at, Phys. 
Rev. Lett. 87, 155001 (2001)]. The EAW is strongest well below the phase-matched frequency for 
electron acoustic scatter, and therefore the EAW is not produced by it. The beating of different 
beam acoustic modes is proposed as the EAW excitation mechanism, and is called beam acoustic 
decay. Supporting evidence for this process, including bispectral analysis, is presented. The linear 
electrostatic modes, found by projecting the numerical distribution function onto a Gauss-Hermite 
basis, include beam acoustic modes (some of which are unstable even without parametric coupling 
to light waves) and a strongly-damped EAW similar to the observed one. This linear EAW results 
from non-Maxwellian features in the electron distribution, rather than nonlinearity due to electron 
trapping. 

PACS numbers: 52.38.Bv, 52.65.Ff, 52.35.Mw, 52.35.Fp 

Keywords: laser-plasma interaction; stimulated Raman scattering; electron acoustic waves; electron acoustic 
scattering; Eulerian Vlasov simulation 



I. INTRODUCTION 

The role of kinetic effects in stimulated Raman 
backscatter (SRBS) is of much current interest. SRBS 
is the three-wave parametric coupling of a pump light 
wave (the laser, labeled mode 0) to a counter-propagating 
daughter light wave (mode 1) and co-propagating elec- 
tron plasma wave (EPW, mode 2), and satisfies the res- 
onance conditions fco = k\ + &2 an( l w o = lo\ + L02- It 
may remove a substantial amount of energy from a high- 
intensity laser propagating through a plasma and create 
energetic electrons SRBS and other laser-plasma in- 
teractions must be sufficiently controlled for laser-driven 
inertial fusion to succeed. In the high-temperature plas- 
mas expected on ignition experiments such as the Na- 
tional Ignition Facility (NIF) 0] and Laser Megajoule 
(LMJ) |3j, linear theory of parametric instabilities pre- 
dicts SRBS to have a weak spatial gain rate due to the 
large feAn (A_d = [t^T e / {nQe 2 )] 1 / 2 is the electron De- 
bye length) and resulting heavy Landau damping. How- 
ever, recent experiments [U, HI, Q and kinetic simulations 
0, H !, M, El El have shown SRBS reflectivities well 
above coupled-mode convective gain levels. 

This kinetic inflation or enhancement of SRBS is 
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widely attributed to electron trapping in the EPW po- 
tential well, which flattens the electron distribution func- 
tion near the EPW phase velocity v v i 0, As has 
long been known, trapping reduces the EPW damping 
rate [3] and induces an amplitude-dependent frequency 
downshift [l5[ . SRBS saturation results from the compe- 
tition of nonlinear effects including the damping reduc- 
tion, frequency shift, pump depletion, and trapped par- 
ticle instability [lGj], and possibly spatio-temporal chaos 
after ion modes build up near the laser entrance [13], 
among others. A theory for the onset of kinetic inflation 
and the resulting time-averaged reflectivity is not yet in 
hand, but would be of great value for designing laser fu- 
sion implosions. 

This current work is aimed at better understanding ki- 
netic effects such as electron trapping in SRBS for time 
scales and frequencies such that ion dynamics can be ig- 
nored. To this end, we have performed 1-D Eulerian 
Vlasov-Maxwell simulations of SRBS using the ELVIS 
code [ill [TH, [r| . A rapid increase of reflectivity above 
the convective value is seen as we vary parameters such as 
the pump intensity. The system does not reach a tempo- 
ral steady state when kinetic enhancement occurs, but in- 
stead the reflectivity comes in sub-picosecond bursts (as 
seen in, e.g., Ref. [19|). Besides enhancing SRBS, trap- 
ping also changes the electrostatic modes present. Elec- 
trostatic spectra from our simulations reveal the EPW 
dispersion curve splits into an upper branch and several 
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beam acoustic modes (BAMs) 0,[T3|, and that SRBS in- 
volves one of the latter. The SRBS daughter BAM activ- 
ity occurs along the streak that phase-matches for pump 
decay, which extends down in frequency from the linear 
phase-matching point with a slope equal to the daugh- 
ter light wave's group velocity v g \ ss — c. This behavior 
has been seen in Thomson scattering spectra from recent 
Trident experiments [2(| and PIC simulations @. 

A low-amplitude acoustic wave (meaning that u> oc fc, 
not an ion acoustic wave), which we call the electron 
acoustic wave (EAW), also appears once inflation devel- 
ops in our simulations. Moreover, weak reflected light 
that phase-matches for scattering off this mode is also 
seen, which we call electron acoustic scatter (EAS). This 
process is reminiscent of what was called stimulated elec- 
tron acoustic scatter (SEAS) and observed in Trident sin- 
gle hot spot experiments 0, Q (although at the most in- 
tense pumps, the experimental SEAS levels relative to 
SRBS are much stronger than in our simulations). To 
our knowledge, this paper is the first to analyze in detail 
EAWs and EAS arising from the SRBS dynamics, and 
not from manually-distorted distributions or other strong 
seeding. SEAS has also been reported in PIC simulations 
of plasmas overdense to SRBS and at relativistic pump 
intensities [2lj . and may have been seen in underdense 
Vlasov simulations [22| . 

Most of the EAW energy in our simulations is much too 
low-frequency to match EAS and must be produced in- 
dependently. We propose the three-wave coupling of two 
BAMs and an EAW, which can be called beam acous- 
tic decay (BAD), as the EAW excitation mechanism. 
Smaller-amplitude EAWs at higher frequency are sub- 
sequently generated by harmonic coupling. The pump 
laser then scatters off these fluctuations, referred to be- 
low as electron acoustic Thomson scatter (EATS). The 
phase-matching needed for BAD and EATS is seen both 
in electrostatic (fc,u>) spectra and bispectral analysis. 

We study the electrostatic mode structure by solv- 
ing the linear dispersion relation for simulation distribu- 
tions projected onto a Gauss-Hermite basis. This yields 
the BAMs and EAW seen in our simulation spectra, 
and shows a damping reduction and frequency downshift 
(compared to the EPW for a Maxwellian) of the BAM 
involved in SRBS. The BAMs are linearly unstable - in 
the absence of parametric coupling to light waves - for 
some wavenumber range. In addition, we find an EAW 
root with heavy damping. This linear EAW differs from 
the nonlinear EAW due to trapping discussed by others 
[23L [24], [25|, [2^, [2?! , which unlike ours is undamped and 
cuts off at k\o ~ 0.53 (in the limit of zero field ampli- 
tude). 

The paper is organized as follows: Section HT1 presents 
the model equations and geometry used in the ELVIS 
code. The dependence of SRBS on pump strength for 
simulations of Trident single hot spot conditions is dis- 
cussed in Sec. IIII1 Kinetically-inflated SRBS, and the 
related beam acoustic modes and electron acoustic scat- 
ter, are carefully studied in Sec. IIVI for a pump inten- 
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FIG. 1: 1-D geometry of the ELVIS code. 

sity of 2 PW/cm 2 . Section fVl demonstrates that similar 
physics occurs in higher temperature, higher density plas- 
mas found in hohlraum fills. We discuss our results (in- 
cluding experimental relevance) and conclude in Sec. IVII 
The Appendix details the Gauss-Hermite mode-finding 
method. 

II. MODEL AND NUMERICAL APPROACH 

We use the Eulerian Vlasov-Maxwell solver ELVIS 
[Til IT2I HH , which resembles the code presented in [28[ . 
Our code implements the simplest model which allows 
for electron kinetic effects in Raman scattering. The ge- 
ometry is one-dimensional, with all spatial variation and 
wavevectors in the x direction. We describe the particle 
species by the nonrelativistic Vlasov equation in x and 
keep the ions fixed in this paper. The light waves are lin- 
early polarized in y, and the particles constitute a cold 
collisionless fluid in this direction (the light waves are 
undamped). The geometry is illustrated in Fig.[TJ 

The governing equations are: 
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m s , Z s e, and n s = J dv x f s are the mass, charge, and 
number density of species s; s — e for electrons, Z e = — 1, 
and e > is the positron charge. A number-conserving 
Krook relaxation operator is included, with relaxation 
rate vrs(x) and equilibrium distribution fo s set to the 
initial Maxwellian (J dv x fo s — 1). We use a large 
Vks = 0.2ti> p (idp = n n e 2 / (e m e )) at the edges of the 
finite length density profile to absorb plasma waves gen- 
erated by SRBS and to reduce fluctuations generated as 
the electron density adjusts to a sheath-like pattern. A 
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nonzero central value of uk s can mimic transverse side- 
loss from a laser speckle, although v Ks = outside the 
damping regions in this paper. 

The Vlasov equation for the distribution function f s is 
solved on a fixed phase-space grid via operator splitting 
[29l ] . The spatial and velocity advections are performed 
via a semi-Lagrangian method (shift along characteris- 
tics) with cubic spline interpolation. The x shift is pe- 
riodic with J s (xl) — fs{xn), although the dynamics in 
the central region are effectively finite in x due to the 
damping regions and sheaths at the edges (subscripts L 
and R denote quantities at the left and right edges, re- 
spectively; see Fig. [1]). For the v x shift, f s is assumed 
zero for \v x \ > «max, which means particles are lost when 
accelerated beyond v ma x- After the v x shift, at each x 
gridpoint we add to f s an initial Maxwcllian with a den- 
sity chosen to cancel the change in number caused by the 
v x shift. The transverse light-wave variables E + and E~~ 
advect to the right and left, respectively. We advance 
E^ without dispersion (in vacuum) by shifts of one x 
gridpoint, which locks the x and t spacings by dx = c dt. 
The extreme edges of the simulation box, beyond the pe- 
riodicity points {xl,xr) of f s , are "moat" regions where 
f s = for all time, but -E propagate until leaving the 
system. We specify E ± (t) at the boundary from which 
each field is advected into the box, and can thereby inject 
a (pump, seed) light wave via (E + ,E~). 

Due to the low numerical noise in our Vlasov code, 
SRBS does not develop if only the pump light wave is in- 
cident and the damping regions are used (although SRBS 
can occur with no damping regions). We therefore in- 
clude a low-amplitude SRBS seed with a frequency that 
satisfies linear phase-matching based on the kinetic EPW 
dispersion relation. The physics is similar in runs where 
the seed was turned off after kinetically-enhanced SRBS 
develops. 



III. INTENSITY SCALING OF SRBS FOR 
TRIDENT CONDITIONS 



Simulations with conditions similar to the Trident sin- 
gle hot spot experiments @ reveal a sudden onset of 
kinetically-enhanced SRBS as the pump strength in- 
creases. A rapid increase in reflectivity was also ob- 
served in the experiments. We use a pump vacuum 
wavelength Ao = 527 nm, background electron density 
no = 0.025n c {uq/uu p — 6.32) where n c is the pump crit- 
ical density (no/n c = (w p /o;o) 2 ), and an electron tem- 
perature T e = 0.5 keV. The central flattop (between the 
Krook damping regions) has length 75.4 /im. We inject 
an SRBS seed with Ai s = 653.4 nm (oji s /oj p — 5.10) 
and intensity Im = 10~ 5 /ol at the right edge, where 
Iol is the pump intensity at the left edge (in this pa- 
per, a constant Iq independent of x and t is under- 
stood to be Iol)- The pump and seed light waves beat 
to drive an EPW with k 2 Xjj = 0.352. In all runs we 
set dx = A2/2O and v x grid spacing dv x < vt r . s /^- 



v tr = 2<x>b//c2 is the trapping island half- width in the 
EPW electric field, ujb = io p {n 2 /no) 1 ^ 2 is the trapped 
electron bounce frequency, and n 2 is the EPW density 
perturbation peak amplitude. Vt r , s is Vt r for the EPW 
produced by beating of the pump and unamplified seed 
light waves; in other words, dv x resolves trapping even if 
there were no seed amplification. The v x grid extends to 
(6.4, 8.6, 9.6)u Te for I 0L (<,=,>)1 PW/cm 2 . 
Fig. [2] shows i? av versus pump strength, where R av 
is the reflectivity R = (Iil) /Iol time-averaged start- 
ing from 1 ps (to allow the lasers to turn on and 
transit the domain) to the end of the run (at 10 ps). 
All runs were below the absolute instability threshold 
Io a = 29.8 PW/cm 2 for undamped light waves, which 
corresponds to the threshold undamped SRBS growth 
rate j 0a = (1/2) /^ g2 1 1 ^ 2 ^2 , with v 2 the EPW Lan- 
dau damping rate. Coupled-mode theory thus predicts 
SRBS approaches a temporal steady state where the seed 
light wave is convectively amplified across the box. The 
steady state can be solved analytically for a three-wave 
model that includes pump depletion but neglects light- 
wave damping, and neglects EPW advection compared 
to damping (the strong damping limit) [30| . The strong 
damping limit is valid for Iql "C Ioa- The scattered light 
wave in this limit satisfies 



dh 

dx 



(7) 



where I = I (x)/I l, h = h{x)/h R , and I - 
(u)q/uji)Ii = constant. L is the length of the gain region, 
and the intensity gain exponent G in practical units is 
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Aip is the scattered light wavelength in the plasma, 
Ao.^i is Ao in microns, and IoL,pw is Iol in PW/cm 2 . 
x{k 2 ,u 2 ) = -(k 2 X D )~ 2 Z'[u> 2 / \k2VTeV2)] is the electron 
susceptibility, where Z is the plasma dispersion function. 
h 2 and ui 2 are given by the beating of the chosen pump 
and scattered light waves. The resulting reflectivity is 



R(l-R + s) = Sexp G(l - R) 



(9) 



with R — (u>o/oji)R and s = (luo/^i)Iir/Iol- 

For weak pumps the simulation gains are small and 
in accord with the coupled-mode calculation, but they 
suddenly rise well above it for Io = 0.4 PW/cm 2 (ki- 
netic inflation or enhancement). i? av barely increases 
with pump strength and saturates around 1% for Io ^ 1 
PW/cm 2 . This is well below the coupled- mode level for 
the largest pump strengths and indicaties nonlincaritics 
besides pump depletion saturate SRBS (kinetic defla- 
tion). 

We do not expect trapping to be important in SRBS 
unless a resonant electron undergoes at least one bounce 
cycle before crossing the gain region. This happens for 
all the runs in Fig. [21 using the convective steady-state 
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I ql (W/cm 2 ) 

FIG. 2: Time-averaged reflectivity i? av and bounce number 
Nb for Trident parameters described in text. Asterisks are 
simulation 7? av . The solid curve is the coupled-mode steady- 
state -R av from Eq. (|9]), and the dashed curve is Nb from 
Eq. (T12"]) ; both include pump depletion. 

fields. The above criterion treats longitudinal endloss 
as the only de-trapping mechanism, since we have not 
included transverse sideloss (e.g., via the Krook opera- 
tor) in our runs. Morales and O'Neil calculate the time- 
dependent damping rate and frequency of an undriven 
EPW as an initial value problem [3, EH- That is, the 
damping and frequency vary with the number of bounce 
orbits the trapped particles have experienced, defined be- 
low as Nb- They find that after the resonant electrons 
complete about one bounce cycle, Landau damping is 
substantially reduced, and the frequency shift is roughly 
its late-time asymptotic value. These results easily trans- 
late to a boundary value problem, appropriate to our fi- 
nite geometry. However, driven EPW's in SRBS may 
differ from Morales and O'Neil's free-wave calculation: 
see the discussion at the end of this section. 

In steady state, resonant electrons emerge from a 
Maxwellian distribution at the left damping region (x = 
0) and undergo Nb{x) bounce orbits (dependent on the 
field amplitude) before reaching the position x. Across a 
domain of interest [0, x], 

N B (x) ee (27T)- 1 f dx' k B (x') (10) 
Jo 

where k B = ^>b/v P 2 is the bounce wavenumbcr in the lo- 
cal field E x (x) . The total amplitude (drive plus response) 
of the EPW driven by two beating light waves as given 
by the linearized Vlasov equation is 

^(W^^. (11) 
no 2 1 + X «j e 

vxe = (Te/m,,) 1 / 2 and v OS i ee eEij > {m e uj i ) is the complex 
electron oscillation velocity in light wave i. A Fourier 



amplitude / is related to its physical field f p by f p — 
(l/2)fexpi(kx~ujt) +cc. In the convective steady state 
described above, n 2 = n 2s (Ioh) 1/2 where ri2s is ri2 for 
= Iol and I\ = I\r. The resulting Nb is 

N B (x) = N Bs {x)~ f dx' (f /i) 1/4 (12) 
x Jo 

where N Bs (x) = [uj p x/{2-Kv P 2)}{n2 a /n ) 1 / 2 is N B (x) if 
there were no seed amplification. Neglecting pump de- 
pletion for illustration, 7j osl — v os i.r cx p[(G/2)(1 — x/L)], 
and Nbr, or the total Nb across the domain of length 
L, is 

e G/4 _ 1 

N BR = Nbs(L) — - . (13) 

Note that Nb s oc depends very weakly on the seed 
intensity. 

Nbr in the convective steady state is shown in Fig. [2] 
as a function of pump strength. Substantial bouncing 
occurs even for the weakest pump (Nbr=1-63 for Io=0.1 
PW/cm 2 ), suggesting that the convective profiles may be 
altered by a damping rate and frequency that vary appre- 
ciably across the plasma with N B {x). However, kinetic 
enhancement only occurs for Nbr ^ 3.5, with A r e,R=3.52 
for the smallest enhanced run (io=0.4 PW/cm 2 ). This 
runs counter to the simple-minded expectation, based on 
the Morales- O'Neil free-wave theory, that SRBS inflation 
should occur for Nbr~ 1 (or even less, since the damp- 
ing rate in Fig. 1 of [15| first becomes zero for Nb ~ 0.6). 
This may indicate a difference between the damping rates 
for free and driven waves. It could also be that inflation 
requires the damping to be small over a sufficiently large 
region, and not just near the right edge. These questions 
merit further research. 



IV. TRIDENT CASE I = 2 PW/CM 2 : KINETIC 
ENHANCEMENT AND ELECTRON ACOUSTIC 
SCATTER 

A. Reflected Light 

We now study in detail the case Iq = 2 PW/cm 2 to ex- 
pose the physics of kinetic enhancement and EAS. Figure 
[3] displays the instantaneous reflectivity, with time aver- 
age R av = 1.03% that greatly exceeds the coupled-mode 
level from Eq. © of 7.75 x 10~ 4 . Moreover, R(t) occurs 
in temporal sub-picosecond bursts instead of becoming 
steady, as has also been reported in PIC simulations (e.g., 
Refs. UM)- 

The time-resolved spectrum of reflected light is dis- 
played in Fig. |3Ja). SRBS is the dominant signal. Early 
in time it occurs at the seed frequency uj\ s but upshifts 
once SRBS becomes strongly enhanced around 3 ps. This 
results from the trapping-induced EPW frequency down- 
shift Scattered light between the SRBS and pump 
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FIG. 3: Instantaneous reflectivity R for I = 2 PW/cm 2 Tri- 
dent run. 



frequencies is observed shortly after the upshift. The 
time-integrated power spectrum, contained in Fig. Hfb), 
reveals EAS is much stronger than the "noise level" in 
neighboring frequencies (the noise is partly due to SRBS 
sidelobes resulting from windowing). The EAS reflec- 
tivity is i? C as=2.06x 10~ 8 , from the reflected light in the 
band w = (5.8 — 6.05)w p during the time EAS is ob- 
served (t > 2 ps). The broad spectrum of reflected light 
for cj < u!i s , larger than the signal above the SRBS peak, 
is scattering off electrostatic fluctuations (possibly from 
BAMs) with lo > L02 visible in Figs.[6](a) and [Til 

All power spectra in this paper were made using 
Welch's method of averaged overlapping modified peri- 
odograms, with Kaiser window functions [U |32| . Ag- 
gressive windowing was needed to reveal weak, time- 
dependent signals, such as EAS in the presence of much 
larger SRBS. The low noise inherent in Vlasov codes al- 
lows such weak processes to be captured. 



(a) Reflected light ( E" ), power decibels (dB) 
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FIG. 4: (Color.) (a) Time-resolved and (b) time-integrated 
spectrum of reflected light E~ for Io = 2 PW/cm 2 Trident 
run. (loo, wis) are the imposed (pump, seed) frequency. The 
lo spectrum of a field / in "power dB" is 201og 10 |/(u)|. The 
EAS asterisk in (b) is at lo = 5.94a;,,. 



B. Electrostatic Activity 

We now turn to the electrostatic activity. Fig. [5] 
presents the envelope of the longitudinal electric field, 
found by rms-averaging E x (x,t) over one SRBS wave- 
length and period. Once SRBS enters the enhanced 
regime, the plasma waves propagate away from the laser 
entrance along fairly well-defined "rays" with group ve- 
locities w vxe- Three broad "bursts," which are strong 
near the laser entrance around 3, 6, and 10 ps, correspond 
to reflectivity bursts in Fig. [3] The initial break-up of the 
envelope into rays, near the laser entrance shortly after 1 
ps, may be due to the trapped particle instability [l6l.l33j|. 
and coincides with the first, isolated spike in reflectiv- 
ity at 1.3 ps. In particular, the envelope is modulated 
at k\o ~ 0.05, which is close to the bounce wavenum- 
ber k B X D = 0.066 calculated for the observed amplitude 
E x e/(m e Lo p VTe) ~ 0.15. Also, the electrostatic spectrum 
over the space-time region of the break-up (not shown) 



possesses sidebands at {k%,u)2) ± (^b,Wb). However, we 
do not see clear signs of the trapped particle instability 
for later times. 

The electrostatic (k,ui) spectrum, taken over all time 
and between the two Krook damping regions, is shown in 
Fig. [51 The "Stokes" curve is the locus of daughter elec- 
trostatic modes (£2,^2) phase-matched for electromag- 
netic decay of the pump, namely, (fej,^) = (ko,&o) — 
(fcijWi), with lui S [Q,Wrj] and ck t — ±(ujf — tUp) 1 ^ 2 
for i = 0, 1. ki > and ki < give, respectively, 
the smaller &2 (forward SRS) and larger &2 (backward 
SRS) legs of the curve. SRBS occurs along the Stokes 
curve, mostly on a frequency-downshifted streak relative 
to the matched EPW. While there is some activity on 
the EPW curve for A; < 0, the curve splits into upper 
and lower branches near k\o — 0.2. The extended activ- 
ity lower in frequency than the EPW curve, and which 
becomes acoustic (u cx k) with a larger slope for small k, 
is analogous to the beam acoustic modes (BAMs) found 
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FIG. 5: (Color.) Envelope (rms) of electrostatic field E x (x, i) 
for 7 = 2 PW/cm 2 Trident run over the x domain between 
the damping regions. 



by L. Yin et al. [9|. 

A low-amplitude acoustic mode {u cx k) with phase 
velocity S3 vx e , which we call the electron acoustic wave 
(EAW), is also present. It is mostly energized for fcArj ~ 
0.15, and remains acoustic but weakens toward higher 
k. The most strongly-excited EAW has a frequency 
~ 0.15o;p well above the ion plasma and ion acoustic wave 
frequencies, so EAW interaction with ions can be ne- 
glected. Spectra over successive time windows reveal the 
EAW is first excited at, and is always strong for, lower k. 
The EAW and EAS both develop after kinetic enhance- 
ment of SRBS begins, which suggests they result from 
SRBS-induccd modifications of the electron distribution. 
The asterisk at (/cAd, w/w p ) = (0.379,0.385) with phase 
velocity 1.02vtc is phase-matched to the EAS asterisk at 
oj = 5.94w p in Fig. [31(b) . However, the EAW curve is 
not generated by EAS, since most of the EAW energy is 
far from the phase-matched point. Also, the amplitude 
is not increased near the phase- matched point, implying 
decay of the pump laser does not significantly change 
the EAW level. The scattering from the EAW appears 
to be off independently-generated fluctuations and can 
be viewed as electron acoustic Thomson scatter (EATS) 
in analogy with Thomson scattering from conventional 
plasma waves. 

We propose that a three- wave interaction involving dif- 
ferent parts of the extended BAM feature generates the 
low-fc EAWs, and dub this process beam acoustic decay 
(BAD). The rest of the EAW curve is weakly excited 
by EAW harmonic coupling (harmonics of an acoustic 
mode lie on the same acoustic curve). BAD is likely 
a two-pump process (instead of a pure parametric de- 
cay), where the BAM curve is energized separately from 
BAD, and BAMs with different (k,oj) then beat to pro- 
duce EAWs. This is supported by the fact that the low-A: 



BAM tail does not have a noticeable peak correspond- 
ing to the daughter BAM in BAD. Moreover, one can 
estimate the ratio of daughter- wave amplitudes expected 
from coupled-mode theory. For weakly damped BAMs 
and a strongly-damped EAW (appropriate for the linear 
modes found in Sec. IIV C|) . and neglecting pump deple- 
tion, the ratio of daughter BAM to daughter EAW ac- 
tions is {ve/lb) 2 where v e is the EAW damping rate and 
7b is the BAD growth rate. For undamped daughters 
the actions are equal (the traditional Manley-Rowe re- 
lation). The E x power ratio from Fig. [5] is very large: 
I E% , b am / E x } e aw | 2 S3 10 4 . Still, without formulae for j b 
and how electric fields relate to actions, it is difficult to 
infer that little of the BAM is produced by BAD. We 
discuss the generation of BAMs below. 

The EAW and resulting EAS occur for moderate pump 
strengths. Neither appear when SRBS is not kinetically 
enhanced. Both increase with Jo above the enhancement 
threshold, until no distinct EAS peak is visible for Iq > 5 
PW/cm 2 . The upshift in uj\ due to the nonlinear down- 
shift of the SRBS EPW grows significnatly with pump 
strength. As the pump increases, there is also more re- 
flected light higher in frcqnuency than SRBS; however, 
this light is spectrally broad and has no peaks (such as 
EAS). A distinct EAW is seen for the I = 5 PW/cm 2 
run but not for Iq = 10 PW/cm 2 . We discuss the re- 
lation of our results to experimental reports of SEAS in 
the conlcusion. 



C. Distribution function and linear modes 

We now examine the electron distribution f e and the 
linear electrostatic modes it supports. Phase-space vor- 
tices develop in the electric field of the SRBS plasmon 
and flatten the space-averaged f e near v P 2- Although 
the vortices are regular for early times, once SRBS be- 
comes enhanced they smear into each other. This indi- 
cates a non-monochromatic spectrum and has been seen 
to result, e.g ., from the trapped particle instability in 
Refs. [ill [Hf. Figure [TJa) displays f e space- and time- 
averaged over periods of 6.3/w p and from 50.6-53.8 fjjn 
(near the domain center) for the Iq = 2 PW/cm 2 Tri- 
dent run. The plateau width is somewhat correlated to 
the local wave amplitude, but f e remains flattened even 
when the field is weak. This suggests nonlocality may be 
important in kinetically-enhanced SRBS. 

The averaged f e , taken from Fig. [7]Ja) at 3.15 ps, is 
plotted in Fig. [TJb) and is fairly representative of f e 
during enhanced SRBS. The observed flattened plateau 
half-width (for the time and space region in Fig. [Hb)) 
of 1.5«T e compares well with w tr = lAlvxe from the ob- 
served amplitude E x e/(rn e ui p VTe) = 0.192. Note that 
f e (vx) in the trapping region is nearly flat (it actually 
increases by 1.6% over its minimum for v x > ^2)- This 
is typical for fully-enhanced SRBS, although in the first 
Raman burst there is a clear period when df e /dv x > 
for v x > v p 2- Such an upturn appears as the pinched 
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(b) E , power dB 




FIG. 6: (Color.) Power spectrum of E x over the (x,t) domain of Fig. [5] for I = 2 PW/cm 2 Trident run. The "EPW" 
curve is the EPW for the initial Maxwellian, and the "Stokes" curve is defined in the text. The "EAS" point at 0.385a;p is 
phase-matched with the EAS asterisk in Fig. |4jb). Panels (b) and (c) are zooms of the SRBS and EAW regions of panel (a). 



bulge in Fig.[7]near (t,v x )—(2 ps, ivre)- Some adiabatic 
calculations of trapping in EPWs show that f e should 
be symmetric about the phase velocity when the wave 
amplitude decreases and electrons are detrapped [3^ |. 

Figure [5] shows Sf = f e — fu at t — 3.15 ps, where 
Jm is a Maxwellian fit to f e for v x £ [— 4,0]uTe (close 
to the initial Maxwellian). Note Sf consists of a broad, 
beam-like structure near v P 2 due to trapping, as well as 
at two lower velocities. Modeling f e as a bulk-and-beam 
distribution, such as a bi-Maxwellian, may not be correct, 
especially for modes like EAWs with low phase velocities 
near the secondary structures. 

To gain insight into the electrostatic dynamics, we 
study the linear modes of the numerically-obtained f e . 
Although a large-amplitude electric field is present and 
may demand a nonlinear calculation, linear theory yields 
good agreement with the observed spectrum. We project 
Sf onto a Gauss-Hermite basis, as described in the Ap- 
pendix. This gives an analytic longitudinal susceptibility 
X, valid in the complex plane, as the sum of x for /m and 
the projected Sf. Figure [5] displays the full Sf and its 



projection up to the N = 12 basis function and with the 
scaled, shifted velocity u = (v/vT e — 3)/0.7. The projec- 
tion only misses fine-scale features which would appear 
in higher-order basis functions. These features should 
weakly distort the real frequency, which stems from col- 
lective oscillations in / e , but may affect the damping rate 
more strongly, which depends on the slope of f e at the 
phase velocity. 

The electrostatic linear modes are obtained by solving 
the dispersion relation e = 1 + x = for complex us given 
real k. We do not include a ponderomotive term due to 
parametric coupling with light waves, which would allow 
for SRBS. Fig. [5] presents the zero contours of e r and e,, 
which intersect at the linear modes e = 0, for f e from 
Fig. [5] This plot over the same (k, us) region only using 
/m bas two zero contours, which intersect at the linear 
EPW; there is an infinity of strongly-damped modes for 
the Maxwellian, which fall below the chosen plot range. 
As N is increased the contours change slightly, and a 
single pair of e r — and = contours can split into 
two pairs. However, the plot remains qualitatively the 
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FIG. 7: (Color online.) f e space- and time-averaged over 
Lu p At = 6.3 and x = 50.6-53.8 Aim for I = 2 PW/cm 2 Trident 
run. Panels (a) and (b) are, respectively, for all times and 
t — 0, 3.15 ps. v P 2 = 3.48uTe is the phase velocity of the 
EPW driven by the beating of the pump and seed light waves. 
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FIG. 8: Sf — f e — fm (solid) and its Gauss-Hermite projec- 
tion (dotted) for (/ e ) from Fig. [7] at t = 3.15 ps. Projection 
parameters given in text. 
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FIG. 9: Zero contours of the electrostatic dispersion relation 
e — for k\o = 0.3, computed from f e as projected in Fig. [8] 



same. This is further discussed in the Appendix. 

Figure [10] graphs five modes, calculated for f e from 
Fig. [8] (t = 3.15 ps), that resemble the observed spec- 
trum. In Fig. [TlJ u r for these modes is superimposed on 
E x (k,u>) computed for t = 2.75 — 4.25 ps; the agreement 
is quite good. The linear modes show the splitting of the 
EPW curve into an upper branch and a set of BAMs, as 
well as the appearance of an EAW. SRBS occurs along 
the mode labeled BAM1, which is frequency-downshifted 
and less damped than the Maxwellian EPW. Thus, much 
of the nonlinear EPW physics of Morales and O'Neil is 
recovered by linear theory with the modified f e , as also 
seen in (35| . 

The upper branch and BAMs are linearly unstable for 
some k, in the absence of coupling to light waves. Since 
the slope of the modified f e is small or slightly positive for 
some v x , it is reasonable to find bump-on-tail or beam- 
plasma instabilities. These beam-driven, growin g el ectro- 
static modes have been studied previously [36l. l37l|. and 
may account for the higher BAM fluctuations compared 
to the level of k < conventional EPWs in Fig. [5] The 
linear instability of BAMs allows them to be excited sep- 
arately from BAD, which based on the discussion at the 
end of Sec. IIVB1 seems to have little effect on the BAMs. 

The linear modes also contain an EAW with u>/k de- 
creasing from 1.05«Te to VTe as k increases. It is always 
heavily damped, unlike the undamped nonlinear EAW 
due to trapping discussed by other workers (e.g., Refs. 
[24l I25I]). An EAW also appears in an analogous mode 
calculation for a bi-Maxwellian f e presented in Ref . [l7| , 
with a beam whose drift velocity is close to that of the 
EAW. It is shown there that stimulated (quasi-mode) 
scattering off the EAW has a small, positive growth 
rate. Possible phase-matched BAD triplets are shown in 
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FIG. 10: Electrostatic linear modes from the projection in 
Fig. |U "EPW" indicates the EPW for the initial Maxwellian. 



Fig. EI Note that the lowest-fe daughter BAMs (white 
circles) lie above the "knee" in the BAM feature, and 
match with daughter EAWs (blue circles) higher in k 
than the strongest EAW activity. The break in the BAM 
streak may thus explain the cutoff in the EAW spectrum, 
although a calculation of the nonlinear coupling mecha- 
nism (e.g., coupling coefficients) would shed more light. 



D. Bispectral Analysis 

To further support the three-wave processes of BAD 
and EATS, we turn to bispectral analysis [38|, |39J. This 
has been used previously to study three- wave interactions 
in p lasma physics, mostly in magnetized plasmas (Ref. 
[40l | is an early example). Although we do not pursue 
it here, bispectral techniques can give estimates of the 
coupling coefficients and mode growth rates [4l|, [52, EH • 

The power spectrum, or the Fourier transform of the 
two-point correlation function, can be generalized to cor- 
relations among more signals at more times, or higher- 
order spectra. Given three real, zero-mean signals x(t), 
y{t), and z(t), their correlation function C 3 and complex 
bispectrum P3 are 

C 3 {x\y,z}(T 1 ,T 2 )= I dtx{t)y(t + n)z(t + T 2 ) (14) 




0C^^ 
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FIG. 11: (Color.) Roots from Fig. [TOl overlaid on E x (k,u) 
for t — 2.75 — 4.25 ps and the x domain of Fig. [6] The white 
and blue circles are possible phase-matched daughter BAMs 
and EAWs, respectively, for beam acoustic decay of a parent 
BAMl root with {k\ D ,u /w P ) = (0.36, 1.1). 



and 



P 3 (lj 1 ,oj 2 ) = / d Tl dT 2 e^^+^Cairun). (15) 



The notation {x\y, z} indicates the fields being corre- 
lated, and Eq. (|16[) shows why x is treated differently We 
compute the bispectrum by a generalization of Welch's 
method for power spectra [32j |. Namely, we divide time 
into Af windows, apply a window function to the signals, 
find their Fourier transforms, and in each window form 

P 3j {u)x,L0 2 ) = Xfai +u) 2 )Y j {u) 1 )Z j (w 2 ). (16) 

P3 = (-P37) where () denotes averaging over the N win- 
dows. P3 measures the signal amplitudes that frequency 
match (as three- wave interactions do). If the signals at 
luq = LO1+U2, lo\, and lo 2 have independent phases <pi but 
constant amplitudes, then the P^^s have random phase 
factors and IP3I ~ A 1 / 2 . This corresponds to three in- 
dependent signals that "gratuitously" frequency match, 
although they are not dynamically coupled. However, 
for a perfectly coherent process like a pure three-wave 
interaction, <fio = 4>\ + 4> 2 , the Pa^'s are in phase, and 
IPs I ~ A. The degree of coherence is quantified by the 
bicohcrence 63: 



(17) 



P 



3jl 



1/2- 



1 63 1 always lies between zero (no phase coherence) and 
one (perfect phase coherence), and its denominator scales 
like N. 

P 3 {E+\E-,E X } is shown in Fig. [U The diagonal 
streak indicates Stokes decay of the pump: lo\ + lu 2 — 
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FIG. 12: (Color online.) Bispectrum \P 3 {E + \E~ , E x }\ for 
t = 5 - 10 ps at x = 52.4 rim for I = 2 PW/cm 2 Trident run. 
The asterisk marks the seeded linear SRBS point, u/i + u>2 — 
ujq along the dominant diagonal streak. 

(a) (b) 
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FIG. 13: (a) bispectrum \Pa\ and (b) bicoherence 1 1>3 1 from 
Fig. 1 121 along the line uii +u>2 = 6.32cj p for decay of the pump 
laser. 



luq = 6.32o-> p . SRBS is clearly dominant, and weaker 
EATS and scattering off BAM noise for uji < uj\ s are 
also present. Figure [TBI displays P3 and 63 along the line 
LU1+LU2 = 6.32cj p ; EATS is seen in both as a pair of peaks 
below 6lu p . Given the incoherent, Thomson- like nature of 
EATS, it is not surprising that it is much less bicohcrent 
than SRBS. 

To examine BAD, we consider the bispectrum of 
three spatial Fourier amplitudes E x (k, i), computed from 
E x (x,t). To distinguish left- from right- moving waves, we 
perform a 2D Fourier transform to the {k,ui) domain, set 
the result to zero for all u> < 0, and then invert the tem- 
poral transform. This procedure can be considered as a 




0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 



°>1 t E x.] 

FIG. 14: (a) Bispectrum and (b) bicoherence of E x (k,t) 
Fourier modes ki\ D = 0.344, 0.209, and 0.136 for i = 0, 1,2, 
over t — 5 — 10 ps in the 7 = 2 PW/cm 2 Trident run. 
E X i = KeE x (ki, t). The peak in IP3I and |&3| shows para- 
metric coupling in the BAD process. 



Fourier transform in x and a Hilbert transform in t. For a 
wave E x = cos(k e x — Lu e t) this gives E x (k e ,t) oc expiu> e i. 
We use E X i = Re E x (ki) to compute P3, which requires 
real fields. 

Figure [T4l presents Pz{E x0 \E x \, E x2 } and the corre- 
sponding £»3, where ko = k\ + k 2 . The time window of 
5-10 ps is when EAW activity is strongest, although sim- 
ilar results obtain for 0-5 ps using slightly different fc's. 
feA_D = 0.13 is chosen slightly below the maximum EAW 
power, and koXo — 0.344 gives a larger P3 and more well- 
defined peak in b 3 than do fc 's closer to the SRBS peak 
of 0.362. The strongly localized peak in P3 indicates 
beam acoustic decay, and persists (although at lower ab- 
solute amplitude) for nearby k choices. A bicoherence 
peak, however, is seen only for certain fc's, which shows 
when the EAW is driven by the beating of the BAMs as 
opposed to satisfying the sum rule with less dynamical 
coupling. The frequencies (ujq, 071,^2) where P3 peaks 
give (k,oj) pairs that lie, respectively, slightly below the 
SRBS peak, near the "knee" in the BAM feature, and 
near the maximum EAW amplitude. 
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FIG. 15: Same as Fig. [2] but for the hohlraum conditions 
described in the text. 



V. HOHLRAUM PARAMETERS 

To demonstrate that kinetically-enhanced SRBS, 
BAMs, EAWs, and EATS are not peculiar to a nar- 
row parameter range, we show in this section that they 
occur in ignition hohlraum fill plasmas as well as the 
single hot spot conditions studied above. We choose a 
pump wavelength Ao = (1054/3) nm, background den- 
sity no = 0.1n c , electron temperature T e = 3 keV, and 
central Krook-free flattop length of 75.2 /jm. The SRBS 
seed light wave has Ai s = 574.8 nm and Im = 10~ 5 io> 
which beats with the pump to produce a plasmon with 
k2^D = 0.357 (similar to the above Trident value). The 
time-averaged reflectivity i? av versus I a is plotted in 
Fig. [T5l and shows an enhancement threshold of Jo = 0.8 
PW/cm 2 , corresponding to Nb — 2.03 bounce orbits in 
the convective steady state. Note that here R av for large 
Jo saturates near 17%, which is much larger than the 
Trident value of 1%. 

We focus on the case Iq — 2 PW/cm 2 . Figure riffla - ) 
displays the time-resolved reflected light spectrum. As 
for the Trident runs, the SRBS light upshifts in fre- 
quency, at which time EAS develops. The SRBS reflec- 
tivity is temporally bursty with large pulses separated 
by about 1.5 ps. The electrostatic envelope (not shown) 
is similar to Fig. [5j with four broad pulses propagating 
away from the laser entrance and correlated to reflectivity 
peaks. The electrostatic spectrum E x (k,uj) is presented 
in Fig. [Tfflb') . Again, the EPW curve splits into an upper 
branch and a set of BAMs, and SRBS is downshifted in 
u> along the Stokes curve. An EAW is present and mostly 
energized at to much lower than the EAS Stokes point. 

As the pump strength I is raised, a distinct EAW and 
EAS persist. This is unlike the Trident runs, where both 
are washed out by very broad EPW frequency downshifts 
for strong-pump cases. In addition, reflected light more 
intense than EAS appears near loq for the hohlraum pa- 
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FIG. 16: (Color.) (a) Time-resolved spectrum of reflected 
light E~ for Iq = 2 PW/cm 2 hohlraum run. u caB — 2.69lu p 
labels the peak of the EAS spectrum, and "re-scat" indicates 
Raman re-scatter of SRBS in its own forward direction, (b) 
Power spectrum of E x over the (a;, t) domain inside the damp- 
ing regions. The stars are Stokes points for EAS with scat- 
tered light of ui/ujp — 2.47 and 2.69. 



rameters (visible in Fig. [TOT a)). and becomes stronger 
with increasing Iq. This cannot be stimulated Bril- 
louin scattering since the ion are fixed, although mo- 
bile ions may significantly distort this signal. The Stokes 
curve passes through a broad, low-frequency feature near 
k\r> ~ 0.45 in the E x spectrum, and the « ojq reflected 
light may be scattering off this. 



VI. DISCUSSION AND CONCLUSIONS 

Kinetic inflation of SRBS and related physics has been 
studied with Eulerian Vlasov simulations in single-hot- 
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spot and hohlraum-fill conditions. Inflation of the reflec- 
tivity above the coupled-mode steady state level develops 
suddenly as the pump becomes stronger, in accord with 
recent Trident experiments. The trapped particle insta- 
bility may modulate and break up the electrostatic field 
envelope early in time but does not dominate the subse- 
quent dynamics. Once trapping flattens the distribution, 
enhanced SRBS occurs in picosecond-scale bursts accom- 
panied by electrostatic pulses propagating away from the 
laser entrance. 

SRBS light upshifts in frequency due to trapping. 
This can be viewed as due to a nonlinear downshift of 
the daughter plasmon frequency, following Morales and 
O'Neil, or as scattering off a lower-frequency linear mode 
supported by the trapping-modified f e (in the spirit of 
quasilinear theory). We adopt the latter approach, and 
find linear BAMs that agree with the observed electro- 
static spectrum and give frequency-downshifted SRBS 
plasma waves. Our Gauss-Hermite linear analysis reveals 
a set of BAMs, some of which are unstable without light- 
wave coupling and thus excitable by purely electrostatic 
dynamics. 

EAWs and EAS have been observed in our simula- 
tions and understood in new ways. The Gauss-Hermite 
method reveals the modified f e supports a heavily- 
damped linear EAW, which differs from the nonlinear, 
undamped EAW reported by Rose and others. It is pre- 
mature to say that EAWs generated by SRBS are a solely 
linear phenomenon or that the nonlinear theories are ir- 
relevant. Linear explanations nonetheless have a certain 
appeal (Occam's razor). Future work should explore the 
impact of de-trapping mechanisms like speckle sideloss, 
and broadband seeding in Vlasov simulations, on EAW 
physics. Preliminary results of simulations we have done 
with a broadband seed show the EAS level is independent 
of bandwidth. This supports our EATS interpretation of 
EAS, as opposed to a stimulated process where seed light 
in that frequency range would be parametrically ampli- 
fied. 

The connection between our EAS picture and exper- 
iment is not yet clear. The relevant Trident experi- 
ments (Ref. Q, Fig. 10) for I = 5 PW/cm 2 saw EAS 
light of 10~ 6 the SRBS intensity, similar to the ra- 
tio in our run with Iq = 2 PW/cm 2 . For pumps above 
5 PW/cm 2 the epxeriments see spectrally-narrow and 
much stronger EAS (~ 1 /3000 the SRBS intensity), while 
our simulations for Iq > 2 PW/cm 2 show no EAS peak 
but instead strongly-upshifted SRBS and broad, higher- 
frequency light (as mentioned at the end of Sec. IIV Bp . 
Moreover, earlier measurements on Trident [4J| with a 
random phase plate (RPP) smoothing the interaction 
beam (giving multiple speckles), which were interpreted 
in Ref. [6| as potentially being SEAS, are much stronger 
than in our runs. Simulations with different density and 
temperature (and thus fc^Au) from this paper, but still 
near the Trident single-hot-spot values, develop EAWs 
and EAS qualitatively like the results presented here (in 
particular, the EAS intensity is ~ 10~ 6 that of SRBS). 



Physics not included in this paper's simulations, such as 
sideloss, ions, or multi-dimensional effects, may provide 
better experimental contact. Measuring the electrostatic 
spectrum in the EAW region of the (k, oj) plane would 
shed light on EAWs and EAS, in particular revealing if 
most EAW activity is at lower fc and u> than the EAS 
matching point (as in our runs). 

Our picture of EAW excitation is beam acoustic de- 
cay (BAD), or the three- wave interaction of two BAMs 
and an EAW. Given the very low EAW amplitudes, this 
is likely a two-pump process, with the daughter BAM 
primarily generated separately from BAD. The daughter 
EAWs can weakly couple to higher-fc harmonics on the 
EAW curve. Moreover, f e is sometimes distorted (but 
not fully flattened) near the EAW phase velocity; this 
may lower its damping rate, or facilitate excitation of 
higher-fc EAWs. The pump laser scatters off the higher-fc 
EAW fluctuations, which we think of as electron acoustic 
Thomson scatter (EATS) . We emphasize that the EAWs 
are strongest well below the matching point for EAS, im- 
plying EAS cannot excite them. In addition, the EAW 
amplitude shows no increase at the EAS matching point, 
which again suggests that EAS is not significantly ener- 
gizing the EAWs. 

After several decades of work, there are many unan- 
swered questions about SRBS. Simple estimates for the 
onset and time-averaged reflectivity due to kinetic infla- 
tion are being investigated but are not yet established, 
and the relative importance of the various saturation 
mechanisms is still unclear. Research on EAWs and EAS, 
besides its intrinsic physical interest, has value in un- 
derstanding SRBS: light scattered off EAWs may indi- 
cate when distribution functions are strongly modified. 
Although large EAS reflectivity has yet to be clearly 
observed, it would not be the first surprise in the rich 
physics of laser-plasma interactions. 
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APPENDIX A: LINEAR MODES VIA 
GAUSS-HERMITE PROJECTION 

The Appendix presents how to find the susceptibility 
and linear modes for an arbitrary distribution by Gauss- 
Hermite projection. We specialize to the 1-D electro- 
static modes of an electron plasma with fixed ions, and 
choose units where u> p =Ai)=VTe = 1 for a reference no 
and T e . The linear susceptibility x is 



X{k,uj) 



dVr, 



dv 



f 



(Al) 



where v p = ui/k is the wave phase velocity, the integral is 
taken along the Landau contour, and / is the background 
distribution with J dv f = n/ng. If we write / = Ylifi 
then x — X)j Xi where Xi is X f° r fi- This is convenient 
for our present work, where / consists of a Maxwellian 
plus a small correction centered at a nonzero v. We use 
a scaled, shifted velocity u = (v — vq)/5v (and similarly 
for Up) and put F(u) = f[v(u)]. We expand F(u) — 

J2 n=0 F n g n (u) over the basis {g n } and obtain x(k,u>) 

and 



■(k 2 Sv) 1 Xu(u p ) where Xu(u p ) 

Xu,n i u p) 



d 
du 



du 



^ FnXu,n('U'p 

g n (u) 



(A2) 



A convenient basis, which gives an inherently localized 
F and analytically known Xu.m is the Gauss- Hermite ba- 
sis 



»(«) 



1 



H n (u)e 



-u 2 /2 



(A3) 



H n is the order n Hermite polynomial (we follow the no- 
tation of [45j). The g n 's are the orthonormal quantum 
harmonic oscillator eigenstates: J_ ao du g n (u)9m(u) — 
S nm . The projection weights F n are given by F n = 
S-oo du 9n(u)F(u). 

The Xu,n' s satisfy a recurrence relation for n > 2, found 
by utilizing properties of H n : 



1/2 



n- 1 



Xu,n{ U p) = _ ( ~ ) Xu,n~l( U p)+ . 

, Th I \ 71 



The base cases n — 0, 1 are 



1/2 



Xu,o(w P ) 



V2 



Z'(w p ); Xu,i 



dXv 



dw 



Xu,n-2{Up). 

(A4) 
(A5) 



w p = Up/y/2, Z'(w) = dZ/dw, and Z is the plasma dis- 
persion function [46(: 



Z{w v 



-1/2 



(A6) 



We evaluate Z with the numerical algorithm in Ref . [47| , 
which provides a high-order rational function approxima- 
tion valid throughout the complex plane. The derivatives 



of Z are a polynomial times Z plus a remainder polyno- 
mial. Therefore, we have 

Xu : n( u p) ~ K z,n+i{up)Z{w p ) + K R:n (u p ) (A7) 



where Kz y n and K Rn are order n polynomials. Eq. (|A4I 
yields a recurrence relation among the polynomials: 
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1/2 



1/2 



Kz.n-l, 
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^Z.n (A8) 



(A9) 



K R , n (u) = -=K z , n - - ^ ;n _ x (A10) 

(All) 
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A" 



fl,n-2 



These formulas hold for n > 2 and allow the polyno- 
mial coefficients to be pre-computed. The base cases are 
K Rfi (u) = -t: 1/4 V2, K RA = -2-k 1 ' a u, Kz i = _ 7r i/4 Mj 
and K z ,2 = ir 1,A ^(l - u 2 ). 

The linear modes are the roots of the dispersion rela- 
tion e = 1 + x(fc,co>) = 0. We find complex u> for real k 
using Newton's method. This requires Xun( u p)i which 
from Eq. (fATf for n > 2 is 



l+2 {u p )Z(wp) + L R . n+1 {u p ). (A12) 



X'u,n( U P ) 



Lz, n and L R ^ n are order n polynomials given by (again 
for'n > 2) 

L z , n +2(u) = -uK z ,n+i(u) + K' Zin+1 (u), (A13) 
-V2K z ,n+i(u)+K'(u). (A14) 

The Gauss-Hermite method has several advantages 
over other approaches. Numerically performing the Lan- 
dau integral for every v p of interest is time-consuming, re- 
quires care in handling the Landau contour, and needs / 
to be analytically continued off the real v p axis. Our tech- 
nique only computes several integrals for the F n 's, auto- 
matically handles the Landau contour via the Z func- 
tion, and is valid throughout the complex plane. Writ- 
ing / as a series of Maxwellians involves some guesswork 
since Maxwellians do not form an orthogonal basis. One 
can use a different scale for H n and the exponential, for 
instance g n (u) <x F[ n (u)e~ u . This gives a simpler re- 
currence relation for Xu,n, but requires more terms to 
reconstruct / well than the g n basis does. We have tried 
both bases and find similar linear modes. 

Unfortunately, the modes given by any Gauss-Hermite 
basis do not converge as N — > oo. For u far inside the 
classical turning points (in analogy with the quantum 
oscillator), g n ~ sin (k„ it -I- <fi n ) where k„ ~ \/2n + 1, 
as a WKB analysis easily shows. Higher-order g n 's thus 
involve oscillation on smaller velocity scales. When ex- 
tended to complex u, these oscillations give rise to expo- 
nential growth. The Xu,n similarly blow up in the lower 
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half- plane (as is well-known for the Z function), with 
|Xun| growing with n for fixed u below the real axis. 
Very small-scale features in / project onto high-n basis 
functions <?„, which produce rapidly growing Xu,n in the 
complex plane and unphysically distort the modes. As 
N increases for the parameters of Fig. the roots vary 



slightly and a single zero-contour of e can bifurcate into 
two. Low-pass filtering before projection could mitigate 
this. Although the Gauss- Hermite method gives much 
insight into the linear dynamics of an arbitrary /, its use 
requires caution and practice. 
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